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1 Introduction 

The Finite Volume method is well adapted to the computation of the solution of pdes 
which are conservation (or balance) laws, for the reason that it respects the property of 
conservation of the pde under study. The mathematical analysis of the application of the 
Finite Volume method to hyperbolic first-order conservation laws can be dated from the 
mid sixties (see jTS62j for example). Concerning the specific problem of the estimate of 
the rate of convergence of the method, the first result is due to Kuztnetsov |Kuz76j . who 
proves that this rate of convergence in L°°{0,t; L^) is of order h^^"^, where h is the size 
of the mesh, provided that the initial data is in BV and that the mesh is a structured 
cartesian grid. Ever since, several studies and results have come to supplement the error 
estimate of Kuznetsov. Before describing them, let us emphasize two points: 

1. The analysis of the speed of convergence of the Finite Volume method is distinct 
from the analysis of the order of the method. In the analysis of the speed of convergence 
of the method, general data (e.g. BV data) are considered. Indeed, here, the problem 
is to show that the Finite Volume method behaves well regarding the approximation of 
the continuous evolution problem in all his features (in particular the creation and the 
transport of discontinuities). On the other hand, in the analysis of the order of the 
method, restrictions on the regularity of the data are of no importance: see the recent 
work of Bouche, Ghidaglia, Pascal BGPO^ on that purpose. 

2. If, numerically, the speed of convergence of the Finite Volume method applied to first- 
order conservation laws is observed to be (at least) of order /i^/^ in the L°°(0, t; L^)-norm, 
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whether the mesh is structured or is not, the preexisting theoretical and rigorous proofs 
of this result appears to be strongly related to the structure of the mesh. 

Indeed, in the case where the mesh is unstructured (what we call a structured mesh is a 
cartesian mesh with identical cells but this can be slightly relaxed |(X;Y98p . the result 
of Kuznetsov has been extended, but to the price of a fall in the order of the error esti- 
mate. Indeed, h^^'^ error estimate in the L°°(0,t; L^)-norm for the Finite Volume method 
applied to hyperbolic conservation laws on unstructured meshes has been proved by Cock- 
burn, Coquel, Lefloch |(X^L94| . Vila jVi]94] and Eymard, Gallouet, Herbin 'EGHOQj for 
the Cauchy Problem {h}/^ error estimate for the Cauchy-Dirichlet Problem [OV04J). We 
emphasize the fact that numerical tests give an order h^/'^ for structured as well as un- 
structured meshes; still, concerning these latter, numerical analysis did not manage to 
give the rigorous proof of the order h^/'^: there is an upper limit at the order h^/^. 

In this paper, we consider the case where the conservation law is linear. More precisely, 
we consider the linear advection problem with a Lipschitz divergence-free speed. 
In the case of if^-data, optimal (with respect to the order of the error estimate) results 
already exist. Under a strict CFL condition, Despres |Des04cj shows an /i^/^-error estimate 
in the Lf^L^-norm for the upwind Finite Volume method applied to the linear advection 
equation with constant speed. (The author deals with the particular case of 2D triangular 
meshes but extensions to higher dimension and general polyhedral meshes are harmless). 
His technique is based on the study of the dissipation of the consistency error by the 
scheme and the adaptation of the Lax Theorem (see also |Des04aj ) . Roughly speaking. 



the consistency error created at time n5t is of order 1 but of order l/V? + 1 after q time 
steps of the scheme. 

On the other hand, we refer to the work of Vila and Villedieu |V V03j , who prove, again 
under a strict CFL condition, an /i^^^-error estimate in the Lf^^ space-time norm for the 
approximation of Friedrichs hyperbolic systems with data by energy estimates (they 
consider explicit in time Finite Volume schemes, for implicit schemes in the case of scalar 
advection equation, see, as they underline it, the result of Johnson and Pitkaranta [JP86] 
who show an /i^/^-error estimate in the L?' space-time norm for data). 

For an initial data uq E BViHf^), we prove in Theorem|21the expected /i^/^-error estimate 
in the L'^L]. norm, under a strict CFL condition in the general case and, under a sharp 
CFL condition if the speed is independent of t. This result is optimal |1'T95| |^ab97| and 
in the context of unstructured meshes and -B\^-data, this is the first optimal result. 
In our proof of Theorem |2l we first show that it is sufficient to consider initial data which 
are characteristic functions. Then we use the same ingredients as Vila and Villedieu 
in |V V03j . The approximate solution satisfies the weak formulation of the problem up to 
an error term corresponding to the consistency error of the scheme. Roughly speaking, 
this error is bounded via the Cauchy-Schwarz inequality by a H^(Ii'^ x [0,t])-like semi- 
norm of the discrete solution denoted Eh{uo,t). A large error means that the scheme is 
very dissipative and consequently, that the approximate solution is smooth, so Eh{uQ,t) 
is small. The /i^/^-estimate follows from the equilibrium between these contradictory 
constraints. 
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The paper is divided into six parts: first we continue this introduction by describing the 
hnear advection problem, the Finite Volume scheme, our results and we give the main 
lines of the proof. In Section |21 we recall some classical results on the Finite Volume 
scheme. In Section El we introduce the weak formulation satisfied by the approximate 
solution. Then we reduce the study to initial data which are characteristic functions and 
we build the test function used in the weak formulation. In Section I5 we prove some 
Energy estimates. The first and second part of Theorem |21 are proved in Sections El El 
respectively. 



Notations 

If {X,fi) is a measurable set with finite (positive) measure and (p G L^{X), we denote the 
mean of over X by 




If X is a set, lx{x) = 1 if x G X, if x ^ X. 

The set of functions with bounded variation in R'^ is the set of functions with bounded 
Radon measures as derivatives: 

BV(R'^) := \u E L\R'^); sup 

where Hv^Hoo := + ' ' ' + Il°°(R'*) fo^' ^ ^ C^(R'^,R'^) and div is the divergence 

operator. The total variation of u G BVlIV^) is given by 

IImIItv := sup 

v5eC|°(R'*,R'*),||ip||oo<i 



u{x)divip{x)dx 



u 



< +00 



u{x)divip{x)dx 



1.1 The linear advection equation 

We consider the linear advection problem in R'^: 

Ut + div{Vu) = 0, X G R"^, t G R+, 

(1.1) 

u{x, 0) = Uo{x), X E fl, 
where we suppose that V G W^'^{'R'^ x R+,R'^) satisfies 

divy(-,t) = VtGR+. (1.2) 

The problem (jl.ip has a solution for Uq G LI^^CR!^); for the purpose of the error estimates, 
we will consider initial data in BV^R!^). 
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Theorem 1. For every Uq G Lj^^(R'^), the problem admits a unique weak solution 

u, in the sense that u G L;^^^(R_|_ x R'^) and: for every G C^(R°' x R+), 

/ / + V ■ V(p)dxdt + / Uo(l){x, 0)dx = 0. (1.3) 

Moreover, we have 

m(x, t) = uoiX{x, t)), V(x, t) G R'^ X R+, (1.4) 

where X G C^(R'^ x R+jR'^) is such that, for every t > 0, X{-,t) : R'^ ^ R'^ is one to 
one and onto. 

LetY{-,t) := X{-,t)-'^, then for every T > 0, X-Id^d andY-Id^a belong to W^'°°{R'^x 
[0,T]). More precisely, there exists Cq > 1 only depending on T, ||F||vfi>°° (^nd d, such 
that, \/x G R'^, Vt G [0,T], 

\\X{x,t)-x\\<Cot, ||VX(x,t)|| <Co, \\dtX{x,t)\\<Co, (1.5) 
\\Y{x,t)-x\\ <Cot, ||Vr(a;,t)|| < Co, \\dtYix,t)\\ < Cq. (1.6) 

Finally, for every t>0, X{-,t) preserves the Lebesgue measure A on IV^, i.e: 

\{X{E,t)) = \{E), for every Borel subset E of K"^. (1.7) 

Proof of Theorem U All the results cited in the Theorem follow from the characteristic 
formula: u{x,t) = UQ{X{x,t)) where X{-,t) = Y{-,t)~^ and Y solves the Cauchy Problem 
dtY{x,t) = V{Y{x,t),t), Y{x,0) = X. To prove the estimates on X, we notice that 
X{x,t) = Z{0;x,t) where Z{T;x,t) denotes the solution of the Cauchy Problem 

— {T;x,t) = V{Z{T;x,t),T), r G [0,t], 
Z(t; X, t) = X. 

We only give the sketch of the proof: global existences for Y and Z and the desired 
estimates follow from the Cauchy-Lipschitz Theorem and the Gronwall Lemma. Besides, 
by ()1.2|1 . the flow preserves the Lebesgue measure on R''. Existence for (jl.lj) follows from 
the characteristic formula and, by an argument of duality, uniqueness also. ■ 

If Uq G L;^q^(R°'), the solution u belongs to C(Il+, Lj^^ili^)). If the derivatives of Uq are 
bounded Radon measures, we have a more precise result. Namely, using the estimates 
above, we obtain 

Corollary 1.1. Let T > 0, there exists a constant Cq>1 depending on \\V\\w^'°°, T and 
d such that, if uq G BV{Yi'^): 

\\u{-,s) - u{-,t)\\Li < Cq\\uo\\tv\s - t\, 0<s,t<T, 
M;t)\\Tv < CoIImoIItv, 0<t<T. 
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The following result is a direct consequence of the conservative property (jl.7p . The case 
f{v) = v'^ will be crucial in the proof of the main result of this paper. 

Corollary 1.2. Let / : R ^ R 6e a measurable function and uq G LI^^CR!^) such that 
f o uq E L^CR^). Then the quantity / /(«(-, t)) is constant. 

1.2 Finite Volume scheme 

The Finite Volume scheme which approximates (jl.lj) is defined on a mesh T which is a 
family of closed connected polygonal subsets with disjoint interiors covering R*^; the time 
half-line is meshed by regular cells of size 5t > 0. We also suppose that the partition T 
satisfies the following properties: the common interface of two control volumes is included 
in an hyperplane of R*^ and 

there exists a > such that I !^r!l^i~J^l'i , ^ i .,r^ ^ (1.8) 

[ \oK\ < a ^h'^ , Va G T, ^ ' 

where h is the size of the mesh: h := s\xp{diam{K), K G T}, \K\ is the d-dimensional 
Lebesgue measure of K and l^i^l is the {d — l)-dimensional Lebesgue measure of OK. If 
K and L are two control volumes having a common edge, we say that L is a neighbor of 
K and denote (quite abusively) L G dK. We also denote K\L the common edge and ihkl 
the unit normal to K\ L pointing outward K. We set 



J K\L JnS 



We denote by Kn := K x [n6t, {n + l)5t) a generic space-time cell, by := T x N the 
space-time mesh and by dK~ the set 

dK-:={LedK,Vl^^<Q}. 

In the same way we set K\Ln := K\L x [n6t, {n + l)6t). 

We will assume that the so called Courant-Friedrich-Levy condition is satisfied: 

^^IVJIlI < (1-0I^I> yKneM, for some e e [0,1). (1.9) 

LedK,- 

Remark 1.1. Under condition (jl.8|) . the CFL condition (jl.9|) holds as soon as 

ll^lloc^t < (l-0«"'^- 

Notice that if ||V^||oo is small, we may choose a large time step 6t. In order to avoid the 
occurrence of terms with factor 6t in our estimates, we add to ()1.9p the following 
condition: there exists Cq > such that 

5t < coh. (1.10) 
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The Finite Volume scheme with exphcit time-discretization is defined by the following set 
of equations: 

u]^ = f Uo{x) dx , MKeT, (1.11) 
Jk 

^ St +W\ E ^^lK-O = 0, VK.gM. (1.12) 

We then denote by Uh the approximate solution of (|1.1|) defined by the Finite Volume 
scheme: 

Uh{x,t)=ul, y{x,t)eKn. (1.13) 



Main results 

From now on, we assume that V e iy^'°°(R'^ X R+, R^) satisfies jOI). We fix a mesh T 
of mesh-size h > satisfying the uniformity conditions ()1.8|1 . We also fix a time step 6t 
such that the CFL conditions (0|) - (irTni) hold. 

We fix two times < t < T and we assume that t = {N + l)6t for some N E N (so 
that Uh{-,t) = '^x&t^k^^^k)- In order to avoid the occurrences of ht~^ terms in the 
estimates we assume that 

h < Cit, for some Ci > 0. (1-14) 

In the sequel, given an initial condition uq G -BV^(R'^), the function u G C(R+, /.^(R*^)) 
denotes the exact solution to p.l|) and Uh its numerical approximation obtained by the 
upwind Finite Volume method (irTT|l - (ITT2|) - (irT!^ . 

In the results and in the proofs, Co > 1 is the constant introduced in Theorem Q and 
Corollary [UH this constant only depends on T, ||l^|| yi/i.oo and d. The letter C denotes 
various constants which are non decreasing functions of a~^, cq, ci, || K||i4/i,oo and d but 
do not depend on h, 6t, ^, uq, t or T. 

Finally, before stating the main results of the paper, we set 

Mn ■.= {KneM : 0<n<N}, 
and we introduce the quantities 



£hiuo,t) := \M;0)\\l2-\M;t)\\l,. 
Remark 1.2. Notice that h~'^Eh{uo,t) is a discrete version of the quantity 



n |2 
K\ 5 



d 



\dtUh{x,s)\'^ + ^ \ Vi{x,s)\\d^Mh{x,s)\'^ dxds. 



i=l 
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We have the following error estimate: 



Theorem 2. Let Uq in BV(R'^). Under a strict CFL condition f M-f^) with ^ > Oj, we 
have : 

\\u{-,t)-uu[-ML^ < CCt'C'\\uo\\Mt'/'h'/' + ej'th), (1.15) 

with = ^. 

Moreover, if V does not depend on the time variable, the estimate is valid with = 1 
uniformly in ^ & [0, 1). 

This result deserves some comments. 

Remark 1.3. If the speed does not depend on x\ V{x,s) := V(s), we have X(x, s) = 
X — jQ\'{r)dr and we can choose Cq = d^^'^ + ||V||oo (which does not depend on T) in 
Theorem^ and Corollarv ll.il 

Remark 1.4. In the proof of the first part of Theorem|21 the dependency on ^ only appears 
in Lemma [5. II 

Remark 1.5 (On the BV-norm of the approximate solution). The possible irregularity of 
the mesh is an obstacle to the existence of uniform (with respect to h) bound on the 
Lji?\4-norm 

E m^L\\ui-ui\ 

of the approximate solution u^- A counter-example to such a result has been given by 
B. Despres in |Des04bj . However, a uniform bound on the weaker norm (notice the weights 

would be enough to prove the error estimate ()1.15|) by use of the techniques of Kuznetsov 
Kuz76j (and, therefore, regardless whether the conservation law under consideration is 
linear, the remark here is true). This uniform bound on Qh{uo, t) is observed in numerical 
experiments and is probably true. However, in this paper, we do not prove this claim and 
it remains an open problem, even in the linear setting. 

Interpolating between the /i^/^-error estimate in the L^L^-norm of |Des04cj and the 
estimate of Theorem |2l we obtain, under a strict CFL condition, a /i^/^-error estimate in 
the L^LP^-noim for data in H^^'P(R^), l<p<2. 

The analogue result for initial data in W^'^ilV^), 2 < p < oo, requires different techniques 
and will be addressed in a subsequent paper. 

1.3 Sketch of the proof of Theorem [2] 

We only sketch the proof of the first part of the Theorem; besides we do not exhibit the 
dependency on t and C, of the error estimates (the constants C may depend on t, T or ^). 
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To prove Theorem |2l we use the fact that the difference between the exact and the 
approximate solutions satisfies the weak formulation f Lemma 13. ip : 

/ {u-uMt + V -^(1))+ [ iu-Uh)i;0)(pix,0)dx = (/i;, + z/;,)(0),(1.16) 

for every G C^(R'^x R_|_), where the error term {fih + i'h){(p) encompasses the consistency 
error of the Finite Volume approximation and depends on the approximate solution. 
We then notice that, by linearity of the equations, and subsequent principle of superpo- 
sition, we may suppose that the initial data Uq G BVlH!^) is the characteristic function 
Uq = 1^ of a set A with finite perimeter. The norm HmoHtv is equal to the perimeter of A. 
In fact, we show that we can also reduce the study to the case where A does not contain 
very thin parts so that, in particular, the volume of the /i^/^-neighborhood of dA 

Ao := {x G R'^ : d{x,dA) < h^^^} . 

is bounded by C\dA\h^^'^. At this stage, we build a function 0o ^ C^iH'^, [—1,1]) such 
that (po = 1 on A\ Aq, 0o = ~1 on A" \ Aq and ||V0o||oo < Chr^l"^ . Then we set 
s) = ^q{X{x, s))l[o,t](s) so that + ^ ■ = on R"^ x [0, t]. After a regularization 
process, the weak formulation yields 

||K-M)(-,t)|Ui < 1^,(0)1 + 

The term /i/i(0) is 

/i,(0) = \K\\u]t'-u],\(l<p-U{.,{n + l)6t)). 

The purpose of the reduction to the case Uq = 1a is an accurate estimation of the term 
/i(0). Indeed, by definition of (j), we have: for < s < t, s) = 1 if x G X{A \ Aq, s) 
and s) = — 1 if X G X{A^ \ Aq, s). More precisely, there exists C > such that, if 

d{Kn, Uo<s<t{X{dA, s) X {s}) > C{h'/^ + h + 5t), 

then is constant on the cell Kn- We denote by M.]^ the set of cells which do not satisfy 
this property; we have 

/x,(0) = Yl \K\\ul^'-ul\(lct^-h{;{n + l)5t)\. 

and the {d + l)-dimensional volume of the cells of A4]^ is bounded by CldAlh^^"^ (Lemma 
13. 5|) . We use the Cauchy-Schwarz inequality and the estimate ||5j0||oo < Ch~^^'^ to get 

<CE{uo,ty/'\dA\y'h'/\ 

A similar technique yields the same bound on Vh{<P)- 
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Thanks to the Energy estimate (Lemma EH)) Eh{uQ,t) < CSh{uo,t), we get 

\\{uH-uK;t)\\L^ < C {\dA\h'/'Sh{uo,ty/' + \dA\'/'h'/') . 

Finally, since the L^-norm of the exact solution is conserved, we get (Lemma 14.11) : 

Sh{Uo,t) < 2\\{Uh-u){-,t)\\Li. 

Together with the preceding inequality, this yields successively £h{uo,t) < CldAlh^^"^ and 
Uu,-u)i;t)U^<C\dA\h'/\ 

2 Classical results 

By continuity in BV of the L^-projection on the space of functions which are constant 
with respect to the mesh T, we have: 

Lemma 2.1. There exists a constant c > only depending on d such that for every 
uo e BV{R'^), 



2 

K£M,LedK 



\uh{-,0) - UoIIli < ca ^ ||mo| 



TV 



Proof of Lemma 12. It Using the weak density of C^(R'^) in BV^R!^), we may suppose 
that Uo e C](R"'). Let K eT, L e dK. Since \x - y\ < 2h for every {x,y) e K x L, we 
have 

\uk-uI\ < / / \uo{x) - uo{y)\dxdy 

2h 



'K JL 

'1 



< - — - — - / / / \Vuo{{l — r)x + ry)\drdxdy. 
\K\ \L\ JkJlJo 

Now we perform the change of variables: (x, y,r) \—>- {w = x — y, z = {1 — t)x + ty,r = r) 
(of Jacobian determinant equal to 1). We have 

\uk-uI\ < 7-^rrjr7 [ \Vuo{z)\( [ [ g{w,z,r)dwdr]dz, 

J B{xK,2h) \Jq JRd / 

where xk G R-'^ is the centroid of K, and g is defined by g{w, z,r) = 1 ii z + rw & K and 
z — (1 — r)w G L, and g{w, z,r) = otherwise. For {z, r) G B{xk, 2h) x [0, 1], we have 
fiid g{u), z,r)dw < 2'^\K\ if r > 1/2 and g{w, z,r)dw < 2'^\L\ if r < 1/2. Finally, we 
obtain 

2'^+^max(|ir|,|L|)/. f 
K-<l < ^T^n / \Vuo{z)\dz. 

I-*^ I l-^l J B{xK,2h) 
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We deduce from the inequalities (|1.8p that for every K ^ T, 



J2 \KiL\\u%-ul\ < 2'^+ia-2 f 

Jb{ 



\VuQ{z)\dz. 



LddK JB{xK,2h) 

Summing on G T, we get 



\Vuo{z)\M2h{z)dz, 



where M2h{z) is the cardinal of the set {K : d{K,z) < 2h}. By p.8|) . we have 

M,,iz)ah'' <j:K:diK,.)<2h\K\ <\Biz,3h)\. 

Thus M2h{z) < ca~^ and we get the first part of the Lemma. 
Let K & T . Similarly, we have 

|u,,(-,0) -Mol < 7-77 / \uq{x) - UQ{y)\dxdy 

K \J^\JkxK 



< / \Vuoiz)\dz. 



'B{xK,h) 

Summing on i^' G T and using the fact that the cardinal of the set {K : d{K, z) < h} is 
bounded by ca~^, we get the second estimate. ■ 

The divergence free assumption ()1.2|) leads to the following identity 
Lemma 2.2. Under the hypothesis on the divergence ofV il.^) . one has 

LedK- L : KedL- 

Monotony 

Under a CFL condition the Finite Volume scheme is order-preserving: 

Proposition 2.1. Under condition the linear application C : (m^^) i— (u]^^) defined 

by M.l^) is order-preserving and stable for the L°°-norm. 



Proof of Proposition EHJ Eq. (fTT^ gi 



ives 



u 



71+1 

K 



LedK- I LedR- 



i.e., under (jl.9j) . u'^^ is a convex combination of u]^, (ul) l^bk- ■ The stability result 
follows from £(1) = 1. ■ 
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L -stability 

From Lemma [2 .21 it is not difficult to see that the quantity Y^xeT conserved. By 

the Crandall- Tartar Lemma |iCT8d] . this fact and Proposition 12.11 imply 

Proposition 2.2. Under condition il.y^) . the scheme C : (u^) ^ (""x"^) ^■^ stable for the 
L^-norm: 



3 Weak Formulation 

In the remainder of the paper, we assume that Uq G BV^R!^), that the function u G 
C(R+, L^{TV^)) is the exact solution to (jl.ip and Uh is the numerical approximation given 
by the scheme (irTT| - (frT^ - (irTT^ . We also assume that (HHl), the CFL conditions (fTlI|l - 
()1.10|) and the condition ()1.14|) are satisfied. 
We intend to prove that Uh satisfies (ll.3|l up to an error term. 

Lemma 3.1. For every (f) G C^(R'^ x [0, +oo)), 



Uh{(j)t + V -^(1))+ Uh{;0)(j){;0) = l^h{(l)) + iyh{<P), (3.1) 

:= Yl mu]t'-um<prK-{^)K{in+i)5t)i (3.2) 

E E K<-ul){VJ^Ml-\K^L\{V-nct^)l^), (3.3) 

where 

JKn Jk JKiL„ 

Proof of Lemma 13. It We develop the first term of the inequality and perform a discrete 
integration by parts: 

/ / Uh{(f)t + V ■V(j))dxdt= E l^l-^^^^/ {(pt + V ■V(f))dxdt 

= E ul l\K\(^{<f))K{{n + l)6t) - {<l))K{n6t)') + ^t\K^L\{V ■ n<j) 

KndM \ LedK 

= J2 (\K\{ul-ulf'){<P)K{{n + l)St)+ J2 m^L\{V-n<prK,u 

K„eM \ LedK 



- J2 \m<{<p)K{o). (3.4) 
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The last sum is equal to J-^^ Uh{-, O)0(-, 0). Rearranging the terms in the second sum and 
using the identity {V ■ n0)^^ + {V ■ = 0, we compute: for every n > 



Plugging this identity and the definition of the scheme ()1.12|) in ()3.4j) . we get ()3.H) . ■ 

We now assume that uq G BVilV^). Before building the test function 0, we show that we 
may reduce the study to more simple initial data. 
Let T' be the cartesian mesh 



k=l 



The diameter of the cells is h' = rf^/^t^/^/i^/^; this mesh is uniform: conditions ()1.8p are 
satisfied for a' = mm(d~'^^'^,d^'^~^^^'^ /2). In particular, by Lemma f2.H we have 

\\vo\\tv < C\\uo\\tv, \\vo - moIIli < C\\uo\\Tvt^^^h^^^, 

where vq is the projection defined by vq := >^ ( / ) These estimates, the L^- 

stability of the scheme and of the evolution equation (jl.ip imply that it is sufficient to 
prove Theorem 121 for the initial data vq. 



Lemma 3.2. \Fed6(A \Bre84 Let vq E BViR"^). For every r/ G R, define 



1 if < ?7 < vo{x), 
Xvo{x,r]) := { -1 if Vo{x)<r]<0, 
in the other cases. 



We have Vq = Xvoi-,v)d'Vj alfnost everywhere. Moreover, 
Jn 



^oIIli = / \\Xvo{-,v)\\L^dri, 



R 



^oIItv = / \\Xvo{-,v)\\Tvdv, 

R 



where the second identity is a consequence of the co-area formula for BV functions. 

Applying this decomposition, we deduce from the linearity of the scheme and of the initial 
problem, that it is sufficient to prove the Theorem for any initial data Uq which is a finite 
sum of characteristic functions 1k'- 

Therefore, from now on, we suppose that uq = 1a with A = Ui<j<iKj and K[,--- ,K[ 
are distinct elements of T' . 
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Let us point out that 

WuqWli = \A\ = lt'^'^h'^'\ 

and since, the boundary OA is a finite union of distinct edges denoted K[\L'i, ■ ■ ■ , K'^\L'^, 
we deduce from the definition of the Total Variation: 

m 

i=l 

With the notations of Theorem ^ we have u = 1b where B C R'^ x R_|_ is defined by 

B := \Jy{A,s)x{s}. 

We are now going to build a test function (f) for (\l.'6\) and (j3.H) . Let T E C°°(R, R) 
such that < r < 1, T{R) = for < 1/3 and T{R) = 1 for > 2/3. We define 
$oeLip(R^R)by 

^ f T{t-^/^h-'/^d{x,dA)) if xeA, 
\ -T{t-^/^h-'/^d{x,dA)) if xeA'. 

Since d{-,dA) is a 1-Lipschitz continuous function, we have ||V$o||oo < ""^^^llT'lloo- 
Let us introduce a mollifier p G C^(IV^, R+) such that J-^^ P = 1 and supp p C -8(0, 1/3), 
we set 00 := t~'^'^h~'^/^ p{t~^/^h^^/'^ ■ ) * "^o- The function 0o belongs to C^(R'^) with the 
bound: 

||V0o||oo < t-''^h-''^\\T'\\^. 

Moreover, we have 

X e A =^ (f)o{x) > 0, and x e A"" =^ (f)o{x) < 0, 

and if we split R"* in three disjoint subsets: R^ = 0o ^({1}) U ^((-1, 1)) U 0o ^({-1}) =: 
A^ U AqU y4_, we have = 1 on A^, = on A_ and 

Ao C {x : d{x,dA) <t'/^h'/^} C UZ^B{xi,l/2f/^h'/^), 



of Aq satisfies: 



where Xi G R is the centroid of K[\L^, 1 < i < m. Consequently the Lebesgue measure 



|Ao| <Cmt^/2/i^/2 = C\\u4Tvt'/^h^l^ (3.5) 

Finally, we define our test-function by 0(x, s) := (pQ^X^x, s))l[o,t](s) for (x, s) G R'' x R_,_, 
where X is the mapping introduced in Theorem ^ Since divV(-,s) = for every s > 0, 
we have 

<pt + V ■V<p =<pt + div(l^0) =0 in R'^ X (0, t). 
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and from the estimate on V0o above and on VX and dfX (|1.5p in Theorem ^ belongs 
to C^(R'^ X [0,t]) and satisfies 

||V0|Uo.(R.,[o,T]) < CCot-'/'h-'/', (3.6) 

\m\L-ii,^.[o,T]) < CC,t-'/'h-^'\ (3.7) 

Moreover, since m(x, s) = uq{X{x, s)), (j) inherits the properties of 0o: for (x, s) G R°'xR+, 

u{x, s) = 1 =^ (f){x, s) > 0, and u{x, s) = =^ (f){x, s) < 0. (3.8) 

Lemma 3.3. We have 



u-Uh){-,t)(f){-,t)- / (M-M^)(-,O)0(-,t) = fih{(i)) + M<P)- (3.9) 

Proof of Lemma 13. 3t The test function (p is not compactly supported and not time 
differentiable at the time t. So it can not be used directly in p.3|l . 1)3.11) . First, remark that, 
since Uq is compactly supported, there exists R> such that suppMft(-, s), suppM(-, s) C 
B{0, R) for every < s < t + 6t, thus may be replaced in p.3|) . ()3.1|) by the compactly 
supported function (p ■ x where x ^ C^(R'^) is such that x = 1 on 5(0, i?). 
To overcome the non differentiability at t we introduce a mollifying sequence to ap- 
proximate the function l[o,t]- Let {4'q)q>i be a sequence of C^(R+) functions satis- 
fying ifjq = 1 on [0,t], ipq = on [t + 6t/q, +oo) and > t/;^ > -2q/6t. We set 
0g(x, s) := (f)o{X{x, s))iljq{s) for every (x, s) G R'' x R+ and q > I. The weak formu- 
lation ()3.1|) reads 

(^j Uh{x,s)(j)o{X{x,s))dx^ ij'g{s)ds + j Uh{-,0)(f)o{-,0) = flh{(pq) + l^h{4>q), 

The sequence of Radon measures {ip'q{s)ds)q converges to —6t, so the first term of the left 
hand side converges towards 

Uh{-,t)(f){-,t) 



R'' 

as q tends to +oo. 

Clearly, for every iC„ G and L G OK', we have {{(l)q)K)q ^ ((^ " n^g)^^)^ 

(l^ ■ n0)^^ and {{(pg) K{n6t))q {(j))K{n5t). Passing to the limit on g, we get 

«/,(-, t)0(-, t) + / Mft(-,O)0o = At/i(0) + '^/i(0)- 
Rd jRd 



Similarly, we obtain from ()1.3p 



u(-,t)0(-,t)+ / u(-,O)0o = 0. 

R'' JR<* 

Subtracting these equalities, we get the result. 

We relate the first term of the equality of Lemma 13.31 to the Lf^L;^,-norm of the error: 
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Lemma 3.4. We have 

||K-n)(-,t)|Ui < / {u-U,,){;m-,t)+C\\Uo\\Tvt'^''h'/\ (3.10) 

Proof of Lemma 13. 4t Recall that from the monotony of the scheme, we have 

< Uh{x, s) < 1, V {x, s) G R"' X R+. (3.11) 

Let us consider the disjoint decomposition R'^ = X{A+, t) UX(ylo, t) UX(A_, t) and notice 
that ()3.8j) . 1)3.111) imply that (m(x, t) — Uh{x, t))(j){x, t) > for every x G R*^. In particular 



< / iu-Uh)i;t)(l){;t). 
Jx{Ao,t) 

Since < u, m/i < 1, we deduce from ()3.5|1 and the conservation property p.7|l that 
/ \{uh-u){-,t)\ < \X{Ao,t)\ = \Ao\ <C\\uo\\Tvt''^h}'^. 

Jx{Ao,t) 

Finally, using the identities 0(-, t) = 1, u{-,t) = 1 on X{A+, t) and 0(-, t) = — 1, u{-, t) = 
on X{A_,t) and the bounds < < 1, we have 



Uhi-,t) - u{-,t)\ = / {u-Uh){-,t)+ {uh-u){-,t) 

X{A+UA-,t) JX(A+,t) Jx(A-,t) 

{u - Uh)i-,t)(j){-,t). 

X(A+uA_,i) 

Summing up this equality and the two preceding estimates, we get ()3.10|) . ■ 
Lemmas 13.31 and 13.41 and the second estimate of Lemma f2. II fand h < Cit) yield 

\\{uH-u){;t)\\Li < l2h{<j)) + M<l>) + C\\uo\\Tvt'/^h'/^. (3.12) 

The first key idea to estimate ^ih{<P) and Vhifp) is to remark that most of the terms in the 
sums of ()3.2|) . ()3.3|) vanish. 

Lemma 3.5. There exists A^jv ^ M.N such that 

f,,{<p) = mup-'-ui)mi-{<p)K{{n+i)st)), 



with the estimate 



5t\K\ < CC^\\u4Tvt'^'^h^'^. (3.13) 
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Proof of Lemma 13. 5t Let us denote by B the open subset of R*^ x [0, t] of all the points 
(x, s) such that s)| < 1. Let Kn G A^at which does not intersect B. In this case, 
is constant (0 = 1 or = —1) on Kn and we have 

m = + m) and VJ^l{<I))1^ = \KL\{V ■ n0)^^, VL G dR-. 

Consequently, the corresponding terms in the definitions of Hh and Uh vanish and the 
equalities of the Lemma hold with 

M], := {KneMN : KnOBy^^}, 

We now estimate the [d + l)-dimensional Lebesgue measure of the union of the cells 
Kn e A^Tv- From the definition of (j), for every Kn E A4]^, there exist (x, s) G Kn such 
that d{X{x,s),dA) < t^/'^h^'^. Since M = Ui<i<mK[ i L^, for every Kn G A<]v> there 
exist {xK,n-, SKr,) G -ft'n and l<i<m such that (i(X(xif„, Si^-^), Xj) < 2t^^'^h^^'^, where Xi 
is the centroid of K[\L[. Thus using the estimates of ()1.6|) on VF and (9^1^, we obtain 

m 

Kn&M]^ j=10<S<t 

C B {Y{xi,n5t),2Cot^''^ h^l"^ + C^dt)) x [n(5t, (n + 

i=l n=0 

and 

m N 

IJ /sTn C U U 5 (y(xi, 2Coti/2 h^/^ + Co5t + /i)) x [n5t, {n + 
Finally, using Cq > 1, 5t < Cq/i and /i < Cit, we get 

4 Energy Estimates 

The two next lemmas contain the second key idea of the paper. The first one shows that 
the L^-error at time t controls £h{uo,t) (the jump of the L^-energy of the approximate 
solution). The second one relates £h{uo,t) to the regularity of the approximate solution. 

Lemma 4.1. We have 

Sh{uo,t) < 2\\{uh-u)i-,t)\\Li. 
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Proof of Lemma I4.lt Using the conservation property of Corollary 11.21 with f{v) = f^, 
we have ||'u(-, t) ||^2 = ||'"o|li2- Thus 

SH{uo,t) = (||w,(-,0)||i. - llwolliO + {\H-Ml^ - \\M;t)\\h) ■ 

Since Uh{-,0) is the L^-projection of Uq on the mesh T, the first term is non positive. 
Finally, < u,Uh < 1 yields u^{-,t) — u\{-, t) < 2 \uh{-, t) — u{-, t)\. m 

In particular, (|3.12|) implies 

< 2/i;,(0) + 2z/;,(0) + C||Mo||Tyt'/V/2. (4.1) 
Lemma 4.2. We have the following identity 



K„&Mn LedK^ \ M£dK- 



1 T/" T/" 

Remark 4.1. We say that the identity of Lemma 14.21 is an Energy estimate, because of 
the analogy with the a priori inequality 

where Uy satisfies the advection-diffusion equation dtu^ + div(V^'Ui^) — z/Amj, = 0. In our 
case the diffusion is due to the scheme (and does not have the regularity and the isotropy 
of the Laplace operator), and the corresponding parameter u is of order h. 

Proof of Lemma 14. 2t Let n > 0. Using Lemma [2.21 to change the index of summation 
in the last sum of the right hand-side of the first equality, we have 



Ei^KO' = E 1^1+ E ^aVH«)'-E E 

Kar KdT \ L&dK- / Ker L-K^dL- 



= E 1^1 + E ^KL^t {uif - E 

I V L(^dK- / L&dK- 

= Ei^i E «^lK)'' 

where we have set for every Kn G M., V^i^K) := dK:^ U {K}, and for every L G dK:^ 



\K\ ' \K\ 
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With this notations, (jl.l2p reads 



LeV~{K) 



for every Kn G Ai, and we have 



El^l E <L(l-<L)Kf- 
■f^6^ [Lev-Cii") L,Mev;;:{K),M^L 



E/ '^KL'^KA4^l'^M 



Using the identity X^/^gv-(ft') a^-i = 1, we compute 

Ei^i(K)'-«T) = Ei^ij E «m( E «^mK-«m) 

= ^El^l E «MmK-«m)'- 
Summing on < n < A^, we get the Lemma. 



5 Proof of the first part of Theorem [21 

In this section, we assume that ^ > 0. We estimate successively Hh{4>) ^^^^ ^h{(f>)- Applyin 
the Cauchy-Schwarz inequahty to the first formula of Lemma f3 .51 we obtain 



m - {(Phdin + my 

St 



On the basis of ()3.7p and the CFL condition ()1.10|) 6t < coh, we have 

\{<P)l - {<P)K{{n + l)6t)\ < \\dt<P\\Lo^6t < CCo6tt-'/^h~'/^ < CCoh^'h-^'^ 
and therefore, by ()3.13p 
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We now write 

+ E E Sii<-^K)\K^L\i{V-i.)lMlL-{V-n<i^)l^), 

=:I + II. (5.1) 
From the Cauchy-Schwarz inequality, we have 

2\ V2 



1/2 / \- \- lT/«JAf/.2 [ 



iii < E,iuo,ty/' [ J2 E \^KLm 



The conditions (jl.8j) imphes 

E \^KL\h <a~^V\\oo\K\ <C\K\, WKnEM. 

Consequently, (|3.6j) and (j3.13j) lead to 

|I| < CCi+^/^E,(no,t)V2||«o|lT''^'/'/^'/'- (5.2) 
We now estimate the term II; notice that it vanishes in case V = Cst. Let Kn G Ai}^ and 

L E dK~, since +- {V ■ n)^^ — V ■ n^L = 0, we have 
Jrl„ 

{V-n)],MlL-{V-n<prR^= I <l>ix,r)i{V-n)l^~V{x,r).nKL)dxdr 

(0(x, r) - {{V ■ n)l^ - Vix, r) ■ ukl) dxdr. 



iRLn 

Using ()3.6p - ()3.7|) and \ {V ■ n)^^ — V{x,r) ■ iiklI < Ch for (x, r) G K\Ln, we get 



Thus, we have 



|II| < CcA 5^ \K^L\hUt\t-'l^h^l\ 

\K„eA1ir \L(idR- ) ) 

As above, for every Kn G we have X^Leex^ < C\K\ and from ()3.13p : 

|II| < CC^-^^\\u4Tvth. (5.3) 
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Summing up (|5.ip - (|5.2j) - (j5.3p . we obtain 

|/i,(0)| + \uHm < C (Co'/'+'i^,K,t)'/lwo||Tyt'/'/^'/' + C^o^iwo||TV't/^) , (5.4) 
and (P?T|) yields 

ShM < C (ci'^^'En{uo,tf''\\uo\\'^yt'l%''' + \\uo\\Tvit"^h^" + Ct+Hh)) . (5.5) 

The following Lemma will allow us to bound Eh{uo, t) by Sh{uo, t). 
Lemma 5.1. 

Proof of Lemma 15. H From p.l2|l and the Cauchy-Schwarz inequality, we have for 
every Kn G A4, 



Due to the CFL condition ()1.9|) . we have J^LeOR- \^kl\^'^ — l-^li consequently, summing 
the preceding inequality on Kn G TWat, we obtain 



|2 



Using again the CFL condition, for every G we have O^ + YliM^dR- ^m^VI-^I) — ^• 
Thus ^Eh{uo,t)/2 is bounded by the first sum of the right hand side of (j4.2|) . ■ 

Plugging this inequality in ()5.5|) . we get 

£,M < c ( cJ/^+^f.(«o,t)n-^/l%llT'?i'/'/^'/' + hohvit'/'h'/' + C^Hh) 



We then use the Young inequality ah < + with 
to derive the estimate 

ShM < cc^+'WuoWMC't'/'h'/' + th). 

Finally, from Lemma I^TD and the inequahties ()5.4|) . ()3.12p . we compute successively 

1^,(0)1 + |z/,(0)| < cc^-^'r'WuoWMt'/'h'/' + e^Hh), 

ii(«,-w)(-,t)iUi < cc^^'r'\\uo\\Mt'/'h'/' + e^Hh). 
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6 Proof of the Second part of Theorem |21 

We suppose that V does not depend on the time variable. If 



sup 



then the proof of the first case allows us to conclude; therefore we can suppose that we 
have \K*\ < 2^^^qj^..- for some K* G T. This cell satisfying JTHI), we 

deduce that the following inverse CFL condition holds: 

h < 2a-^\\V\\ooSt. (6.1) 

We now bound fih{(p) and z//i(0). From Lemma and the estimate ()3.7|) on dt4>, we have 

Since V does not depend on t, we have = V^, for every Kn G ^A. Consequently the 
discrete time derivative := m^) satisfies the scheme fll.l2|l . The L^-stability 

of the scheme and Lemma 12.11 imply 

t s-^ ,,,,, 1 n , . . ...Q IL.O „,0 



J2 \K\\ul^'-ul\ < <t$^|\^, 
K„eMN KeT KeT 



KLW'^K "^L 



< t\\y\\ooJ2\KiL\\u''K-ul\ <C\\uo\\Tvt. 
KeT 

Thus (using the CFL condition 6t < Coh), we have 

M)\ < C\\uo\\Tvt'/'h'/'. (6.2) 
As in the proof of the first case (see (jSH])), we write z^h(0) = 1 + 11. Recall that 

and that, from the estimate (|3.6|) on V0, we have 



\l\<CCoh'/H-'/' E E ^^I^klIKl-u 



K\ 



Using the definition of the scheme ()1.12|) . for Kn ^ M.\j and L G dK^ , we have 

\K\ 

MedK~ 
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Therefore, 



2|„,n 



L,MedK~ 



Using again the L^-stabihty of the scheme, we have 

|Ii| < CllwollTy 
For I2, we use the Cauchy-Schwarz inequahty to get 

/ 



II2 



< 



\Kn&MN L,MedK- 



T/n T/" f)f^ 

\K\ 



1/2 



M 



U, 



1/2 



E E 

K„&Mn L,MedK~ 



yn yn 

\K\ 



1/2 



1/2 



Thanks to Lemma 14 .21 the first term is bounded by \^£h{uo, 

and by ^J^, the 

second term is bounded by CCQ^'^\\uo\\Tyt^^'^h^^'^St~'^^'^. Finally, using the inverse CFL 
condition II2I is bounded by CCQ^'^\\uo\\Tvt^^'^^'^^'^^h{uo,ty^'^ and we have 



(6.3) 

The proof of the estimate ()5.3|1 of the previous Section is still valid and we have |II| < 
uoWrvth. Summing up this estimate, ()f).2j) and (jfi.Hj) . we obtain 

\^h{(p)\ + 

< CC'J'^' [\\uo\\Tv{t'/'h'/' + cS/'th)+S,{uo,ty^^^^^ . (6.4) 

Plugging ()6.4|) in ()4.H) and using Young inequality to absorb the term £h{uo, t) in the left 
hand side, we are led to 

£h{uo,t) < CC^-^'Wuohvit'^'h'/^ + th). 

The second part of the Theorem follows from the last estimate, ()6.4p and ()3.12|) . ■ 
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